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ABSTRACT 

In the context of the brane-world scenarios with compactified extra dimensions, we 
study the production of brane fluctuations (branons) in hadron colliders (pp, pp and e ± p) 
in terms of the brane tension parameter /, the branon mass M and the number of branons 
N. From the absence of monojets events at HERA and Tevatron (run I), we set bounds 
on these parameters and we also study how such bounds could be improved at Tevatron 
(run II) and the future LHC. The single photon channel is also analyzed for the two last 
colliders. 



PACS: 11.25Mj, ll.lOLm, 11.15Ex 



1 Introduction 



Since rigid objects do not exist in relativistic theories, it is clear that brane fluctuations 
must play a role in the so called brane world, proposed some years ago by Arkani-Hamed, 
Dimopoulos and Dvali (ADD scenarios [1]), where the Standard Model (SM) particles 
are confined to live in the world brane and only gravitons are free to move along the 
D > 4 dimensional bulk space (see [2] for recent reviews). This fact turns out to be 
particularly true when the brane tension scale / (r = f A being the brane tension) is much 
smaller than the D dimensional or fundamental gravitational scale Mp, i.e., f « Mp. 
In this case the only relevant low-energy modes of the ADD scenarios are the SM particles 
and branons which are the quantized brane oscillations. Branons can be understood as 
the (pseudo)Goldstone bosons corresponding to the spontaneous breaking of translational 
invariance in the bulk space produced by the presence of the brane. It has been pointed 
out that branons could solve some of the problems of the brane-world scenarios such as 
the divergent virtual contributions from the Kaluza-Klein tower at the tree level or non- 
unitarity of the graviton production cross-sections [3] . As Goldstone bosons, branons are 
in principle massless, but in the cases where the metric of the extra dimensions cannot 
be factorized, they can become massive [4, 5]. This is similar to the case of pions which, 
being the Goldstone bosons of the spontaneous breaking of chiral symmetry, acquire some 
mass due to the explicit breaking of the symmetry induced by the quark masses. 

In previous works the different effective actions have been obtained, namely: the 
effective action for the SM fields on the brane, that for the branon self-interactions and 
finally that corresponding to the interaction between SM fields and branons [4]. In general, 
this branon effective action can be parameterized by the number of branons N, the tension 
scale / and the branon masses (for an explicit construction see [6]). Using the effective 
action it is possible to obtain the different Feynman rules, the amplitudes and finally the 
cross-sections for branon production from SM particles. In [7, 8, 9] the case of electron- 
positron colliders has been considered. By using the Large Electron-Positron Collider 
(LEP) data it is possible to set important bounds on the tension scale and on the branon 
mass for a given branon number. Other restrictions have also been set from astrophysical 
and cosmological considerations due to the fact that branon dark matter can present 
relevant abundances [10]. 

In this work we study branon production in hadron colliders and also in electron-proton 
colliders such as HERA. Most of these cross sections have been studied by Creminelli and 
Strumia for the massless branon case [9]. We reproduce their results and extend the 
analysis for an arbitrary branon mass. The paper is organized as follows: In Sec. II we 
shortly review the branon effective action. In Sec. Ill we consider the case of proton- 
(anti)proton colliders like Tevatron or the future Large Hadron Collider (LHC). In Sec. IV 
electron(positron)-proton colliders like HERA are studied. In Sec.V we show the main 
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results for the relevant examples and in Sec. VI we set the conclusions. 



2 Effective action 

The relevant effective action describing the low-energy interactions of SM particles and 
branons was derived in [7] , where the necessary vertices are detailed. The branon effective 
action can be expanded according to the number of branon fields appearing in each term: 

S eff [K] = Sfj f [K]+S%[K] + ... (1) 

where the zeroth-order term is just a constant and the second-order is just the free action: 

SfjM = \j d A x{5 aP d^d^ - JW>V), (2) 

with n a (x) the branon fields where a — 1,2, N and is the squared mass matrix 
which, without loss of generality, can be assumed to be diagonal. The effective action for 
the SM particles and their interactions with branons is given by 

Ssm, = J d A x[C SM + ^(4^^°^/ - V ^M 2 a ^ a ^)T^ M ], (3) 

where Csm is the SM Lagrangian and Tg V M is the SM energy-momentum tensor defined 
as: 

n v M = -brtsM + 2^)u=^, (4) 

where is some arbitrary metric on the world brane and rj^ is the Minkowski metric. 

In this work we are interested in the interactions between quarks and gluons or photons. 
Thus, for Dirac fermions with masses m q belonging to some representation of a gauge 
group, such as U{\) em or SU(3) C , with generators T a , the Lagrangian is 

C q = q(iYD^ ~ m q )q, (5) 

where the covariant derivative is defined as = <9 M — hA^T a , h being the appropriate 
gauge coupling. Thus the energy-momentum tensor is given by 

= l -(q(l»D u + YD^)q - (D v qY + D"qY)q) 

-r)^(i(qYD p -D p qY)q-2m q qq), (6) 

from where it is possible to find vertices such as irnqq and irnqqA. For gauge fields A the 
appropriate Lagrangian for perturbation theory is: 

£a = -\f^FZ u + £fv, (7) 



where as usual F® v = d VL A a v — d v A a il — hC ahc J^' il A c v and C^-p is the Fadeev- Popov Lagrangian 
including the gauge fixing and the ghost terms. The energy-momentum tensor is: 

n u = f^f^w + \ v p wn + Tpp, (8) 

from where we can obtain the 7T7tAA, 7nrAAA and ttttAAAA vertices. 

Therefore, by using these energy-momentum tensors and the effective action above, 
it is possible to obtain the different Feynman rules involving branons. One important 
observation is that in all the vertices obtained above, branons appear always by pairs. In 
fact they interact in a way similar to gravitons since they couple to the energy momentum 
tensor. This can be seen by making the formal identification of the graviton field h^ v which 
appear in linearized gravity with 

V - ~\{^d^ a d v ^ - ^M^ttV 3 ). (9) 
«J 4 

where k = Ayfir/Mp and Mp the Planck mass. Of course the physical meaning is com- 
pletely different for branons and gravitons. In any case branons are expected to be weakly 
interacting and then they will scape to detection. Hence their typical signature will be 
missing energy and momentum. Since branons are produced by pairs, the energy spectrum 
of any other particle present in the final state will be continuous. In the following sections 
we will study the production mechanisms relevant for the different hadronic colliders. 



3 Proton- (anti)proton colliders 

For the case of proton-antiproton colliders like Tevatron, the most important processes 
for branon production are quark-antiquark annihilation or gluon fusion giving a gluon 
and a branon pair; and (anti)quark-gluon interaction giving an (anti) quark and a branon 
pair. Therefore the expected experimental signal will be in both cases one monojet J 
and missing energy and momentum. This is a very clear signature that in principle 
can be easily identified. Another potentially interesting process is the quark-antiquark 
annihilation giving a photon and a branon pair. In this case the signature is one single 
photon and missing energy and momentum. 

The Feynman diagrams contributing to the main subprocesses qq — > gim, gg — > g^n, 
qg — > qirn and qg — > qmr are shown in Fig. 1, Fig. 2 and Fig. 3. 

From these diagrams and the Feynman rules coming from the effective action of the 
previous section, it is possible to obtain the differential cross section: 
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Figure 1: Feynman diagrams associated to the qq — > g7T7r subprocess. 



where s = (pi + P2) 2 , * = (Pi — <?) 2 ) u = (p 2 — q) 2 and A; 2 = (&i + k 2 ) 2 , P\ and p 2 being 
the anti-quark and quark four-momenta respectively, q the gluon four-momentum and k% 
and hi the branon four-momenta. We have assumed for the sake of simplicity that all the 
branons are degenerated with a common mass M and that all the quarks are massless. 
We have also neglected the effects of the top quark. In addition we have the well-known 
relation s + t + u = k 2 . The contribution to the total cross section of the process pp — > grni 
coming from this subprocess is given by 



where qp(y; s) and q p (x; s) are the distribution functions of the anti-quark inside the an- 
tiproton and of the quark of flavor q inside the proton respectively, and x and y are the 
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fractions of the proton and antiproton energy carried by the subprocess quark and anti- 
quark. The different limits of the integrals can be written in terms of the cuts used to 
define the total cross-section. For example, in order to be able to detect clearly the mono- 
jet one must impose a minimal value for its transverse energy Et and a pseudorapidity 
range given by rj min and r] max . Then we have the limits k min = AM 2 , k 2 max = s(1—2Et/ Vs) 
and t min ( max ) = -(s - k 2 )[l + tanh {r] min ( max ))}/2. On the other hand x min = s min / s and 
Umin = x m i n /x where s is the total center of mass energy squared of the process and 

s min = 2E 2 + AM 2 + 2E T ^/E 2 + AM 2 . (12) 

In addition the dots in (11) represent the contribution of the case in which the quark 
comes from the antiproton and the anti-quark comes from the proton. 
The cross-section of the subprocess gg — > gnu is given by 
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Figure 3: Feynman diagrams associated to the qg — > qirir subprocess. The qg —> qirir subprocess has 
the same diagrams, but changing the quark lines by the corresponding antiquark ones. 



dcr{gg -> gnu) _ a s N(k 2 - AM 2 ) 2 I AM 2 

dk 2 dt ~ 40960/% 2 s 3 te V*" 
(s 4 + t 4 + u*-k s + 6k\s 2 + t 2 + u 2 ) - Ak 2 {s 3 + t 3 + u 3 )), (13) 

where the Mandelstan variables are defined as in the previous case, with pi and P2 being 
the initial gluon four-momenta, q the final gluon four-momentum and k = k\ + k% the 
total branon four-momentum. Then the contribution to the total cross section from the 
PP ~ > 5' 7r7r reaction is 

Vgg(pp -> gmr) = j dx f dyg(y; s)g(x; s) 

f-dk 2 dt d ^2 9nn) ■ (14) 

Jk 2 mm Jt min dk 2 dt 
Here g(x; s) is the gluon distribution function of the (anti)proton, x and y are the fractions 
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of the proton and antiproton energy carried by the initial gluons and the integration limits 
remain the same. From the above equations, it is possible to compute the total cross- 
section a(pp — > gnir) in terms of the cut in the gluon (monojet) transverse energy E T . 
For the qg — > qirn process the cross-section is given by 



da{qg -» qnn) _ a s N jk 2 — AM 2 ) 2 _ UP 

dk 2 dt ~ 2 184320/ W*uV k 2 [uk +Ats ^ 2uk + s h (15) 
with pi and p 2 being the quark and the gluon four-momenta respectively, q the final 
state quark four-momentum and k\ and k 2 the branon four-momenta. The Mandelstam 
variables are defined as in previous cases. The cross-section for the conjugate process 
qg — > qTTTT is exactly the same. Then the total cross section for the the reaction pp — > qmr 
is 

a(pp -> qTTTT) = j dx j dy ^g{y\ s)q p (x; s) 

k2 -dk 2 dt M^2^ + ... (i6) 

J twin CtrC (JAj 

mm iivi.il. 

In this equation x and y are the fractions of the proton and antiproton energy carried 
by the subprocess quark and gluon. The different integration limits are defined as in the 
previous case in terms of the minimal transverse energy of the quark (monojet) Ex and 
the dots refer to the case where the initial gluon is coming from the proton and the quark 
is coming from the antiproton. In addition we have the contribution from the conjugate 
case where we take an anti-quark from the proton and a gluon from the antiproton and 
conversely. This amount just to a factor of two. 

From all the above equations it is possible to compute the total cross-section a{pp — > 
Jtttt) in terms of the cut in the jet transverse energy Ex- 

For the subprocess qq — > ^tttt the cross-section is given by 

Here the notation is similar to the qq — > gmr case with the obvious differences in couplings, 
color and charge factors. Thus 

a (jpp _^ 77r7r ) = j dx j dy^ qp(y, s)q P (x; s) 



f-dk 2 dt**^ 7 ^ + ... 

ib 2 if '->->■ 



min 



dk 2 dt 

All the previous discussion about branon production in pp reactions can be easily trans- 
lated to the pp case. The only point is to change the antiproton distribution functions of 
the different partons by the corresponding proton ones. 
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Figure 4: Feynman diagrams associated to the qq — ► jirir subprocess. 

4 Electron(positron)-proton colliders 

For electron(positron)-proton colliders like HERA, the most interesting branon creating 
process is branon photoproduction, where a photon emitted by the electron(positron) 
interacts with a quark(antiquark) from the proton giving a quark (antiquark) and a branon 
pair. Thus the experimental signature is again one monojet J plus missing energy and 
momentum. The relevant Feynman diagrams are shown in Fig. 5 and the corresponding 
differential cross-section for the subprocess 'jq — > qim is 



da(<yq -> qnn) 3Q 2 q aN{k 2 - AM 



r2\2 



dk 2 dt 



184320 f 8 n 2 sHu 



AM 2 



{uk 2 + Ats 2 )(2uk 2 + t 2 + s 2 ), (19) 



where s — (jp + q) 2 ,t = (jp— k) 2 and u = (q — k) 2 , p being the photon, q the proton quark, 
q' the final quark and k the total branon momenta respectively. The total cross-section 
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Figure 5: Feynman diagrams associated to the qj — > qmr subprocess. The 57 — > g7T7r subprocess has 
the same diagrams but changing the quark lines by the corresponding antiquark ones. 

for the process e^p —>■ qirn is given by 



a(e ± p 



qmr) 



I dx f dy^F(y)q p (x]§) 



k 2 

"max 



dk A 



Umin 
imax 



dt 



da(jq — > qmr) 
dk 2 dt 



(20) 



x and y are defined in this case as q = xP p and P = yP e with P p and P e being the proton 
and electron(positron) momenta respectively. Thus at high energies compared with the 
proton mass s = xys where s = (P e + P p ) 2 - The integral limits y m in, x min , k 2 ninmax and 
t m in,max are defined like in the proton- (anti)proton collider case. 

The photon spectrum F(y) can be obtained from the well-known Weizscicker- Williams 
approximation [11]: 



a 
2ny 



[l + (l-2/) 2 ]log 



Ami 



(21) 



with s' = xs and m e being the electron mass. 
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The cross-section a(e ± p — > g7T7r) can be obtained in a similar way. Then the total 
contribution to monojet plus missing energy and momentum production for large enough 
Et coming from branons can be written as the sum of a(e ± p — > qmr) and a(e ± p — > g7T7r). 



5 Results 

By using the cross-sections shown in the previous sections it is possible to compute the 
expected number of branon pairs produced in the different hadron colliders in terms of 
the brane tension parameter /, the branon mass M and the number of branons N. To 
this end we have used the distribution functions which can be found in [12]. The values 
of the electromagnetic and strong couplings have been taken at the electroweak boson 
masses, namely a = 0.0781 and a s = 0.1171. However our final results do not depend 
too much on the precise value of these couplings. In fact our main source of error is the 
use of an effective action to describe the SM particles and branon interactions since, in 
principle, this is only guaranteed for energies well below 4.7r/. 
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Figure 6: Exclusion region at the 95 % c.l. in the parameter space / — M for TV = 1 from ZEUS data 
on jet production. 



As discussed in the introduction, our main goal in this work is to study the bounds 
that can be set on the /, M and N parameters coming from hadron colliders. We will 
present all our limits at the 95% confidence level. In particular, for the electron(positron)- 
proton case, HERA is the most relevant experiment. In fact, the ZEUS collaboration has 
studied the jet production in charged current deep inelastic e + p scattering. Its results 
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are perfectly compatible with the SM background and therefore, we can set some bounds 
on the branon production and hence on the /, M and N parameters. These data were 
taken from 1995 to 2000 at a maximum CM energy of 318 GeV. The total integrated 
luminosity was 110.5 pb _1 and the cuts on the pseudorapidity and the transverse energy 
were — 1 < rj < 2 and Et > 14 GeV (see [13] for more details). By using the same cuts 
with our cross- sect ions for monojet plus a branon pair production, we find the bound 
/ > 16V 1 / 8 GeV for massless branons. For a branon mass larger than 152 GeV there is 
no restrictions on the / value because of kinematical reasons. For the intermediate M 
values the bounds obtained can be seen in Fig. 6 where we have assumed N = 1. For 
other N values one just has to take into account that the bound scales like V -1 / 8 since 
all the cross-sections are proportional to Nf~ 8 . 
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Figure 7: Exclusion region / — M at the 95 % c.l. for N = 1 with D% data in the monojet channel, and 
with CDF data in the single photon channel. The dashed line corresponds to the LEP-II limits obtained 
by the L3 collaboration using single-photon data [8] 

In the pp case the most relevant experimental information so far is the one obtained 
at the Tevatron (Run I). The -D0 detector has studied the monojet channel and CDF 
the single photon one. As far as the number of events found in both cases is compatible 
with the SM background, we can set new bounds on the branon theory parameters. For 
light branons the most important bound comes from the D$ data. These data were 
taken from 1994 to 1996 at a CM energy of 1.8 TeV and correspond to an integrated 
luminosity 78.8 ± 3.9 pb~ . The cuts on the pseudorapidity and the transverse energy 
were \r)\ < 1 and Et > 150 GeV (see [14] for the details of the analysis). The total 
number of monojets observed was 38 and the expected number from the SM plus cosmic 
rays events was 38 ± 9.6. By using our cross sections for monojet plus a branon pair 
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production with these cuts we get the bound / > 157 TV 1 / 8 GeV for light branons. The 
restrictions for / improve up to a branon mass of 822 GeV. For the intermediate M values 
the bounds obtained can be seen in Fig. 7 for N = 1. 




200 400 600 800 1000 

M (GeV) 

Figure 8: Sensitivity estimation at the 95 % c.l. for the second run of Tevatron in the parameter space 
/ - M for N = 1. 

In a similar way we can use the CDF data on single photon production. In this case 
the total luminosity collected was 87 ± 4 pb -1 and the pseudorapidity cut was \q\ < 1. 
For the transverse photon energy several cuts were considered (for example 55 GeV at the 
75% efficiency). The total expected background for this process was 11.0 ± 2.2, without 
taking into account the QCD contribution (see [15] for the details of the analysis), and the 
number of events found was 11. Comparing this result with our computations for photon 
plus one branon pair production, we find the bound / > 148iV 1//8 GeV for massless 
branons and no bound for M larger than 872 GeV. The bound obtained for the rest of 
the cases is shown also in Fig. 7. 

In addition to this analysis corresponding to the Tevatron data (Run I), it is also 
interesting to make some estimation about the bounds that could be set from future 
experiments such as Tevatron (Run II) and the LHC. In the case of the Tevatron (Run 
II), which is already in progress, the main novelties are a CM energy which equals 1.96 
TeV and an expected integrated luminosity Cji at the end of the run of about 200 pb -1 . 
The detectors are also improved so that the pseudorapidity cuts can be taken as \q\ < 3 
for D® and \rj\ < 3.6 for CDF. This would result in a factor of \J Cn/ Cx on the statistical 
significance when compared to the Run I, with integrated luminosity Cj, provided that 
the CM energy and the cuts were the same. For massless branons, the bound on / scales 
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Figure 9: Sensitivity estimation at the 95 % c.l. for the LHC in the / - M plane for N = 1. 
as the CM energy Eq^. Even more important is the possibility of exploring higher branon 



masses, since the kinematical limit is given by M = yE CM — 2E C mE t /2. In Fig. 8 we 
show the expected bounds from the Run II in the f — M plane, again for N = 1. 

The LHC will produce pp collisions at a CM energy of 14 TeV and the integrated 
luminosity will be something about 10 5 pb _1 . In order to estimate the bounds on the /, 
M and N parameters that will be possible to obtain at the LHC, we have proceeded in a 
similar way as in the Tevatron case, with the obvious changes in the distribution functions 
due to the fact that now we are dealing with pp instead of pp collisions. We have kept the 
same cuts except for the transverse energy which has been corrected in order to maintain 
the same proportion relative to the CM energy. Again the best bounds for / come from 
monojet production, which for M = turns out to be / > lOTSiV 1 / 8 GeV. For low / 
the best bound for M is given by the single photon channel (Mo = 6781). The LHC 
sensitivity for other values in the f — M plane can be found in Fig. 9 for N = 1. 

6 Conclusions 

In this work we have studied the flexible brane-world scenario, where the brane tension 
scale / is much smaller than the fundamental D-dimensional gravitational scale Mb- In 
this case, the relevant low-energy degrees of freedom are the SM particles and the brane 
fluctuations or branons. From the corresponding effective action, we have calculated 
the relevant cross-sections for different branon searches in hadronic colliders. We have 
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used the information coming from HERA and the first Tevatron Run in order to get 
different exclusion plots on the branon mass M and the tension scale / plane for a given 
branon number N . Monojet production turns out to be the most efficient process for light 
branons, whereas the single photon channel is the most important one for heavy branons. 
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87 


55 


-1.1, 1.1 


1.3 10- 10 


148 
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Teva-II 1 


1.96 
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Teva-II 2 


1.96 


10 3 


55 


-3.6, 3.6 
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LHC 1 


14 


10 5 
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-3, 3 


1.8 lO" 11 


1075 


6481 


LHC 2 


14 


10 5 


430 


-3.6, 3.6 


3.8 10~ 12 
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6781 



Table 1: Summary of the main characteristics of the analysis performed for hadronic colliders. All the 
results are presented at the 95 % c.l. We have studied two channels: the one marked with an upper 
index 1 is related to monojet production, whereas the single photon is labelled with an upper index 2 . 
We considered four different experiments: HERA, the I and II Tevatron runs and the LHC. Obviously 
the data corresponding to the two last experiments are estimations, whereas the first two analysis have 
been performed with real data, y/s is the center of mass energy associated to the total process; C is 
the total integrated luminosity; Et is the transverse energy cut; rj m i nymax , the pseudorapidity limits; <jq 
is the estimation for the cross section sensitivity limit; /o, the bound in the brane tension scale for one 
massless branon (N = 1) and M the limit on the branon mass for / = 0. 

We have also extended the analysis to future hadronic colliders. The corresponding 
sensitivity regions for the second Tevatron run and the LHC have also been plotted (see 
Table 1 for a summary of the analysis). 

These analysis improve those already done for electron-positron colliders for heavy 
branons, whereas for light branons, the results are similar [7, 8, 9]. The Tevatron (run I) 
limit M = 872 GeV can be compared to the analogous limit from LEP II M = 103 GeV 
[8]. According to the previous estimations, the Tevatron run II could also improve the 
bound /o = 180 GeV obtained by LEP-II. On the other hand, LHC could detect branons 
up to a mass of several TeV (M = 6781 GeV) improving even the CLIC prospects 
(M ~ 2500 GeV) [7]. 

The study of branons in colliders can be complemented with other bounds coming 
from astrophysics and cosmology (see Fig. 10). In fact, as shown in [10], the branon 
relic abundance can have cosmological consequences. Other issues related to branon 
phenomenology, such as their radiative corrections to the SM processes, or their distinctive 
signatures at colliders with respect to the KK gravitons will be analyzed elsewhere. 
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Figure 10: Relic abundance in the f — M plane for a model with one branon of mass: M. The two 
lines on the left correspond to the fls r h 2 — 0.0076 and flsrh 2 — 0.129 — 0.095 curves for hot-warm relics, 
whereas the right line corresponds to the latter limits for cold relics (see [10] for details). The lower area 
is excluded by single-photon processes at LEP-II [8] together with monojet signal at Tevatron-I. The 
astrophysical constraints are less restrictive and they mainly come from supernova cooling by branon 
emission [10]. 



Note added: After this paper was completed, we were informed by M. Spiropulu 
that CDF collaboration had performed a monojet study [16] which could improve the 
bounds in a more detailed analysis. 
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